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Avalanches of rearrangements in quasi-2D emulsion hopper flow 

Xia Hong/' Kenneth W. Desmond, b Dandan Chen c ' d and Eric R. Weeks* 17 


We experimentally study the flow of a quasi-two-dimensional emulsion through a constricting hopper shape. Our area fractions 
range from <f> = 0.83 — 0.99, such that the droplets are always in contact with one another and are in many cases highly deformed. 
At the lowest flow rates, the droplets exit the hopper via intermittent avalanches. At the highest flow rates, the droplets exit 
continuously. The transition between these two types of behaviors is a fairly smooth function of the mean strain rate. The 
avalanches are characterized by a power law distribution of the time interval between droplets exiting the hopper, with long 
intervals between the avalanches. There is little or no dependence of the flow behavior on the area fraction of the system. 


1 Introduction 

Many slowly strained materials exhibit intermittent flow be¬ 
havior: long still periods punctuated by rapid avalanches 
where material flows—' —. Examples include diverse phenom¬ 
ena such as earthquakes^ 7 , general deformations of solids—, 
stick-slip friction due to granular layers^S, Barkhausen noise 
in magnetic materials—, and sheep herded through constric¬ 
tions^.. More commonly, avalanches are seen in slow flows 
of athermal soft materials such as emulsions^, bubble rafts—, 
foams-i 14-17 , and granular materials^ - -? 3 . These soft materi¬ 
als typically have amorphous structure, necessitating that flow 
and rearrangements are disordered on a microscopic scale. 
The slow flow speed is a key feature: for example, a rotating 
drum experiment with sand inside demonstrated avalanches at 
low rotation rates and smooth flow at high rotation rates?!. For 
granular materials, static friction can prevent the material from 
flowing and can lead to avalanches. In systems composed of 
fluids such as foams and emulsions, stresses are supported not 
by static friction but rather surface tension, which resists the 
deformation of the bubbles or droplets. 

In this manuscript, we present an experimental study of a 
flowing emulsion where we see a range of flow behaviors. 
Our emulsions are oil droplets in water and are compressed 
between two parallel glass plates so that the droplets are de¬ 
formed into pancake-like disks. The area fractions are all 
above jamming—, such that the droplets touch each other and 
are in many cases highly deformed by their neighbors. At the 
slowest flow rates, we see large avalanches of rearrangements. 
At higher flow rates rearrangement events are always occur¬ 
ring and droplets exit continuously. Intriguingly, the transition 
between the two flow behaviors occurs fairly smoothly as the 
flow rate is increased, and at moderate flow rates we see an 
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intermediate type of flow behavior. Despite the range of area 
fractions we study (<|> = 0.83 — 0.99), the area fraction does 
not seem to control the flow behavior. 

2 Experimental Method 

2.1 Design 

Our emulsions are mineral oil droplets in water using Fairy 
detergent (mass fraction 0.025) as a surfactant to prevent co¬ 
alescence of the droplets 25 i 26 . The droplets are produced us¬ 
ing a standard co-flow micro-fluidic technique-? 7 . The radius 
polydispersity of our droplets made this way is 1% (standard 
deviation divided by mean). To prevent droplets from orga¬ 
nizing into crystalline arrays, for each experiment we make a 
bidisperse emulsion by mixing together two separate batches 
of monodisperse droplets at a volume ratio of about 1:1. While 
each individual batch of monodisperse droplets has a low 
polydispersity, there is some variability between batches. The 
mean diameter of the large droplets is 270 ± 50 /./m and of 
the small droplets is 200 ± 40 /an, and the diameter ratios 
of the bidisperse mixtures we form are in the range d^/ds = 
1.5 ±0.2. 

In our experiment, we confine droplets between two 25 mm 
x 75 mm glass slides. The slides are separated by pieces of 
100 /rm transparency film sealed with epoxy. These pieces of 
film act as spacers and thus creating a gap between the slides. 
This gap ranges from 115 to 140 fim in different experiments. 
Nonetheless, within a given sample chamber, this gap is con¬ 
stant with uncertainty 1.8% within any given experiment so 
the slides are parallel (the corresponding maximum angle be¬ 
tween two slides is 1°). The diameters of the oil droplets are 
chosen to be larger than the gap of the sample chamber. Thus, 
the droplets are squeezed between the two glass slides without 
overlapping to achieve a quasi-2D system. 

The left panel in Fig.[j]shows the schema of the chamber. 
The pieces of film are cut to form a symmetric hopper channel 
with angle 0 = 54 ± 5° (see Fig. QJ and opening width 0.7 — 1 


1-011 






Fig. 1 Schema of the sample chamber (left) and raw image of the 
emulsion (right). 


mm. The sample chamber is tilted at an angle 5 ± 1° relative 
to the horizontal, to use the buoyant force of the droplets to 
balance the viscous friction between droplets and glass slides. 
The buoyant force is due to the density difference between wa¬ 
ter and mineral oil (pwater = 1-00 g/cm 3 , p Q ;i = 0.83 g/cm 3 ). 
First we load the emulsion into the sample chamber, and then 
behind the emulsion we add pure mineral oil. A syringe pump 
injects additional mineral oil into the chamber at constant flux 
rate to push the emulsion through the chamber and thus fun¬ 
nel the droplets through the hopper exit. We use a microscope 
with a 1.6x objective lens to image the system, focusing on 
the chamber midplane where the 2D droplet images are clear¬ 
est. A CCD camera records the images in the region close 
to (0.5-2 mm away from) the hopper opening. Depending on 
the mean speed of the flow in a given experiment, the camera 
frame rate is between 0.2 and 2 images/second. This is suf¬ 
ficient to track the trajectory of each individual droplet using 
standard software-^, even at the maximum velocity 0.06 (D)/s, 
where ( D ) is the mean diameter of the droplets. The right 
panel in Fig.Q]shows a typical raw image, in which we record 
hundreds of droplets within the field of view. 


2.2 Control Parameters 

The two control parameters for our experiments are the area 
fraction (|) and the flux rate F. (|) is the fraction of the area 
occupied by oil droplets, which is calculated based on the 2D 
image. We take a total of 45 data sets with 0.83 < (|) < 0.99 
and 0.0001 < F < 0.02 ml/hr. From the post-processed im¬ 
ages, we observe that (j) has only minimal fluctuations during 
an experiment, with a relative standard deviation no more than 
0.5%. These fluctuations are primarily due to the finite field of 
view, with (|) changing when droplets flow in and out. In flow¬ 
ing suspensions of solid particles there can be a self-filtration 
effect—, but in our data the mean area fraction as a function 
of time has no monotonic ascending and thus there is no self- 
filtration. There is an additional possible systematic uncer¬ 
tainty for 4> as the apparent size of each droplet depends on the 
illumination settings of the microscope. We keep these set¬ 
tings constant between each experiment. For each experiment. 


F is set by a syringe pump and thus is constant. However, the 
observed flow velocity fluctuates. This is likely due to some 
compliance in the sample chamber, allowing sample to flow in 
slightly without having to flow out, and building up pressure 
until it is released by droplets flowing out. 

To simplify the discussion, rather than focusing on flux, we 
use the experimentally measured quantity of the mean strain 
rate (y). (y) is defined and calculated as follow, x and y are de¬ 
fined using the coordinate system on the right panel in Fig. [T] 
For each droplet, we determine the trajectory of its center of 
mass, [^(f),y(f)]. Using a short time interval, we also deter¬ 
mine its instantaneous velocity (v x ,v y ). For a given droplet, 
we consider the motion of it and its neighbors (defined as 
those droplets connected by a Delaunay triangulation of all 
the droplets). Using these data, the instantaneous strain rate of 
this droplet j is averaged over all neighbors: 


Vj = l( 


V x (j)~ v x( 0\ 

/! ! 

yj-yt 


(i) 


where the subscript i indicates the ith neighbor of the reference 
droplet j, and the average is taken over all of the neighbors. 
The strain rate can be interpreted as the spacial gradient in y 
direction of the velocity in x direction; the gradients of the y 
velocity are smaller and accordingly we neglect them. The 
mean strain rate (y) of each data set is calculated by averaging 
over all droplets and all time. 


3 Results and Discussion 

We observe a wide range of flow behaviors as we vary F and 
(]) for different experiments. For large F, droplets flow con¬ 
tinuously and smoothly (referred as smooth flow cases in this 
article). For small F, we see avalanche-like flow (referred as 
avalanche cases). For intermediate flux rates F, we observe 
intermediate cases between these two flow patterns. As will 
be discussed below, we do not see any clear dependence of 
these flow patterns on the area fraction <j). 

We summarize these three flow behaviors in Fig. [2] The 
three pictures in Fig.|2ja)-(c) use color to show the time each 
droplet exits the hopper opening to the right. Red droplets exit 
the earliest, and blue the latest. The left picture is a smooth 
flow case, which shows a smooth gradient in color. The right 
one shows an avalanche case, where droplets have distinct 
groups of colors indicating that droplets exit the hopper in 
bursts. Note that the color scale of each plot corresponds to 
a different amount of time, as specified in the caption. 

Figs.[2jd)-(f) quantify these pictures by showing the cumu¬ 
lative number of droplets that have exited the hopper as a func¬ 
tion of time for our three flow cases. In the smooth flow case 
(d), the data form a smooth curve with a well-defined slope, 
showing that droplets exit the hopper continuously at a fairly 
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Fig. 2 Description of the three flow behaviors, (a-c) Images of the samples at a particular time, with the color indicating the time when the 
droplet exits. The red droplets exit earlier and blue droplets exit later. The time ranges for the colors are (a) 1250 s, (b) 3000 s, and (c) 

14000 s. (d-f) The number of droplets that have exited the hopper as a function of time, (g-ij Histograms of the number of droplets exiting the 
hopper within a short time window, chosen such that the mean of the histogram is 10 droplets. The flow conditions are: (a,d,g) smooth flow, 

0 = 0.87, (y) = 0.065 s _I . (b,e,h) Intermediate, tj) = 0.96, (y) = 0.02 s _I . (c,f,i) Avalanche, <j) = 0.96, (y) = 0.01 s _l . 
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Fig. 3 Schema of the definition of interval At. The left figure is at 
time t\ when the black droplet exits the hopper. The right figure is at 
time ?2 when the red droplet exits the hopper. 


constant rate. The intermediate case (e) shows fluctuations 
in the rate, although it is still fairly continuous. In avalanche 
case (f), there are stretches of time where no droplets exit, 
followed by discrete sudden flow events where many droplets 
exit within a short period of time, indicated by the vertical 
portions of the data in (f). Specifically, the first vertical line at 
t « 6000 s relates to all of the yellow droplets in (c) that exit 
at nearly the same time. Again, the existence of avalanches 
despite the constant flux set by the syringe pump shows that 
there is some compliance in the plumbing, such that the pres¬ 
sure builds up before an avalanche. 

Fig. |2g)-(i) show the histograms of exiting flux. The ex¬ 
iting flux size is the number of droplets exiting during a time 
window T, with T chosen to make the mean flux size to be 
10. The smooth flow case (g) has a Gaussian shape while the 
avalanche case (i) has a few rare but large events. To quan¬ 
tify this, the skewness values for these distributions are (g) 
0.15, (h) -0.03, and fi) 2.2 for smooth flow, intermediate, and 
avalanche cases respectively. Not surprisingly, the avalanche 
case has a large positive skewness, and this is generally true 
that all avalanche flow cases have positively skewed distribu¬ 
tions. Given that the avalanche cases have very few events 
overall, our skewness data are noisy and we cannot resolve 
any clear trend in the skewness as a function of our control 
parameters. The general picture shown in Fig.[3g)-(i) is clear, 
though, that avalanche cases have distributions with positive 
skewness and there is a trend toward more symmetric distri¬ 
butions with skewness ss 0 as F increases. 

To better quantify the difference of these flow behaviors, 
we introduce a method focusing on the temporal behavior of 
the flow. In avalanche cases, discrete sudden flow events are 
separated by time intervals where droplets barely move and 
no droplets exit the hopper. Accordingly, we define the time 
between two successive droplets exiting the hopper as the in¬ 
terval At. As shown in Fig. [3] we set t\ as the time when the 
black droplet exits the hopper, to as the time when the next 
droplet (in red) exits, and then At = to — 1\. 
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Fig. 4 Typical examples of three types of probability distribution 
functions of At: (a) exponential distribution, (b) intermediate case, 
(c) power law distribution. The area fraction c) and mean strain rate 
(y) are as indicated in each panel. In (a) the line shows an 
exponential fit P(Af)~e~ A, / T with x = 6.3 s. In (b) the straight line 
is a power-law fit P(A?)~Af~ a with a = 2.1 and the curved line is 
an exponential fit with x = 12.7 s. In (c) the line is a power-law fit 
with a = 1.6. 
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It is apparent in the plots in Fig. [2fd-f) that the distribu¬ 
tions of At are different for the smooth flow and avalanche 
cases. In smooth flow, the values of At are small and do not 
fluctuate much. In the avalanche case. At is sometimes small 
(vertical portions, where many droplets exit over a short time 
interval) and sometimes large (horizontal stretches, where a 
long time passes between one droplet exiting and the next). 
Figure |4]shows the probability distribution functions for At for 
the same three data sets shown in Fig. [2] The smooth flow case 
shown in Fig.[4}a) is well fit to an exponential, as shown by the 
dashed red line; note this is a semilog plot. The exponential 
fit P(Ar)~e~ A, / T defines the mean time between events x, and 
the fit suggests that the time between events follows a Poisson 
process, where events occur continuously and independently 
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Fig. 5 Phase diagram of fitting patterns of P(At) in terms of area 
fraction <) and mean strain rate (y). Red circle: power law; blue 
triangle: intermediate; black cross: exponential. 
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with a constant mean rate. The avalanche case shown in panel 
(c) is well fit to a power law, as shown by the dashed red line; 
note this is a log-log plot. The fit in this case is given by 
P(Af)~Af -ct with « = 1.6, and the power law regime covers 
more than 2 decades in At and more than 4 decades in proba¬ 
bility. The tails correspond to the long periods of time where 
droplets barely move. The intermediate case in panel (b) is 
plotted on log-log axes, and can be fit with either a power law 
(straight line) or an exponential (curved line); neither fit is per¬ 
fect. The exponential fit fails for the largest At while the power 
law is not adequate to describe the small At region. 

In our experiments we have varied both (|) and flux rate 
(quantified by the mean strain rate (y)). For each experiment, 
we use the shape of P(At) to describe its flow behavior. Fig¬ 
ure [3 shows the phase diagram of fitting patterns. There is 
no obvious trend with (]), but more clearly a transition from 
avalanche flow (red circles) to avalanche flow (black cross) 
with increasing (y). Note that the judgment about the best fit¬ 
ting function is done by eye. The quality of each fit depends on 
which range of data is used for the fit, and while we have tried 
several ways to approach the fitting procedure more systemat¬ 
ically, none seem satisfactory for the intermediate cases, and 
none affect the appearance of Fig.0in any substantial way. 

The phase diagram of Fig. 0 is perhaps unsatisfying as the 
intermediate cases (blue triangles) are mixed in with the other 
two cases. However, by ignoring <\> and focusing only on the 
flow rate dependence, the data become more unified. In par¬ 
ticular, Fig. 0a) shows the relation between the power law 
exponent of P(At) and (y). The power exponent a increases 
as the mean strain rate increases. Even when the power law 
fit is not perfect (triangles), the data still follow the general 
trend started by the well-fit power law cases (circles). Smaller 
values of a indicate a broader distribution, where the large At 


events are more significant: these are the avalanche cases with 
long pauses between short bursts when many droplets exit. 
This is similar to previous experimental studies of sheared 
granular materials, where they have power law distributions 
of various stick-slip event properties including forces, energy, 
and avalanche sizes 21 - 30-34 . 

To consider the exponential cases, we note that the expo¬ 
nential fitting parameter X corresponds to the mean interval. 
Accordingly, in Fig. 0b) we compute (At) for every experi¬ 
ment (no matter the distribution shape) and see that this de¬ 
creases as mean strain rate (y) increases. This is plausible 
since the strain rate is bigger for higher flow velocity and thus 
droplets exit faster, which leads to small mean interval. Note 
that for those cases in Fig. 0a) with a < 2, it is problematic 
to define (At) since the mean of a power law distribution with 
a < 2 is infinite. Thus, our finite values of (At) merely reflect 
the finite amount of data. It is plausible that the power law 
distributions may well have some cutoff at very large At. To 
an extent, the relationship shown in Fig. 0b) is trivial by di¬ 
mensional analysis: (At) and (y) 1 have dimensions of time. 
The dashed red line in Fig. 0b) is a power law fit with ex¬ 
ponent — 1.2, fairly close to the -1 suggested by dimensional 
analysis. The different exponent, and the reason the data don’t 
completely collapse, is that the two quantities are different av¬ 
erages: in general (x) ^ (x -1 ) -1 and in specific (At) focuses 
on droplets exiting whereas (y) focuses on the behavior of all 
of the droplets everywhere, averaged over both time and space. 

4 Conclusions 

We have demonstrated that in a single experimental geometry, 
we see behaviors changing from clear avalanches to smooth 
continuous flows as we increase the mean flow rate by a fac¬ 
tor of 100. We quantify these behaviors by examining the 
distributions of times At between subsequent droplets exiting 
the hopper. Intriguingly, the transition in the flow behaviors 
is somewhat smooth as we increase the flow rate: the power 
law exponent characterizing the tails of P(At ) smoothly varies 
as the flow rate increases past the point where a power law 
no longer adequately describes the data. One possibility is 
that at any flow rate, the distribution P(At) may be describ- 
able by a power law with an exponential cutoff, and this cutoff 
may smoothly move to smaller At as the flow rate increases. 
However, the data we have for the intermediate cases [such as 
shown in Fig. 0b)] are hard to interpret in the tails, and so 
it is difficult to resolve this question. The rate dependence of 
our observations is consistent with prior studies of athermally 
sheared 2D amorphous solids which demonstrated rate depen¬ 
denc e- 1 ^ — —. The dependence on velocity is also displayed in 
experimental works on sheared granular materials, where fric¬ 
tion plays a key role^IS. For hopper flow in granular exper¬ 
iments, the presence of static friction can make jamming and 
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Fig. 6 (a) The power law exponent a as a function of mean strain 
rate (y). For the power law fits, only data in the tail are used for the 
fit. Different choices of the minimum At used for the fit give rise to 
different values of a, reflected in the error bars shown here, (b) 
Mean interval (At) as a function of (y). The dashed line is a fit with 
slope of -1.2. 


a T1 event 42 -. So, for fast flow rates in our experiment, a 
rearrangement event near the hopper exit which allows one 
droplet to exit may not be fully completed before the next re¬ 
arrangement happens, allowing for substantial cross-talk be¬ 
tween given events^. At slower flow rates, individual events 
can be less correlated, allowing fewer droplets to exit at a time, 
until the stress builds up and is released in a large avalanche 
event where many rearrangements happen nearly simultane¬ 
ously. In our intermediate cases, it may be that in some mo¬ 
ments rearrangement events occur close enough in space and 
time to influence and enhance each other, whereas in other 
moments events are more individual; some evidence of this 
has been seen in simulations of 2D foam flow— and other sys¬ 
tems^. 

In summary, we see that the flow of an emulsion through a 
hopper can vary from avalanche-like to continuous. The tran¬ 
sition between these behaviors is not abrupt, but rather a con¬ 
tinuous function of the flow rate. At the lowest flow rates, the 
power law exponent we observe approaches a = 1, showing 
that the flow has extremely long quiescent intervals in between 
the avalanches. The decrease of the power law exponent with 
decreasing flow rate [Fig. [6f a)] suggests that even with these 
slow flows, we are not in a quasi-static limit, in agreement 
with a prior study of slowly sheared bubble rafts-ii. In this 
simple limit where the strain rate approaches zero, the flow 
is not simple, but rather dominated by the rare intermittent 
avalanches. 


clogging obvious, where stress-supporting solid arches form 
across the exit 41 . In addition to static friction, such experi¬ 
ments are also driven by a constant force (gravity), whereas in 
our experiments the syringe pump increases the pressure until 
flow occurs, and so no arches can persist indefinitely. 

The transition from avalanche to smooth flow happens 
around (y) ps 0.02 s , indicated in Fig. [5] While we use 
mean strain rate to indicate each data set, actually the strain 
rate is spatially inhomogeneous. Within the imaged field of 
view, y near the hopper exit can be 5 to 14 times bigger than 
the mean strain rate (y). Thus at the hopper exit using this 
local y, the transition point is y« 0.1 s _1 to 0.28 s 1 . The 
inverse of this gives us a time scale Xy = 1/y about 3.6 to 10 
seconds. Xy is the time for the system to strain by 1, which 
microscopically can be viewed as the time for one droplet to 
slide past a neighboring droplet. In a 2D system, the simplest 
topological rearrangement is the neighbor exchange of a group 
of four droplets, known as a T1 event 17 ~ 6 : two droplets that 
are neighbors move apart, and between them two droplets that 
weren’t neighbors before move together and become neigh¬ 
bors. The time scale of a T1 event Xri in our system is a few 
seconds^. This is determined from the mean structural re¬ 
laxation time based on the change of local stress field around 
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